Self-consistent Fermi surface renormalization of two coupled Luttinger liquids 
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Using functional renormalization group methods, we present a self-consistent calculation of the 
true Fermi momenta k p (antibonding band) and k b F (bonding band) of two spinless interacting 
metallic chains coupled by small interchain hopping. In the regime where the system is a Luttinger 
liquid, we find that A = k b F — kp is self-consistently determined by A = Ai[l + g'o ln(Ao/A) 2 ] -1 , 
where go is the dimensionless interchain backscattering interaction, Ai is the Hartree-Fock result 
for kp — kp, and Ao 2> A is an ultraviolet cutoff. For <?o m (Ao/Ai) 2 3> 1 even weak interachain 
backscattering leads to a strong reduction of the distance between the Fermi momenta. 
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I. INTRODUCTION 

What survives of the Luttinger liquid properties of a 
one-dimensional metal if two or more metallic chains are 
coupled by some weak interchain hopping t±? Motivated 
by the proposal that the anomalous normal state proper- 
ties of the high-temperature superconductors are a man- 
ifestation of a two-dimensional non-Fermi liquid statei, 
many authors tried to find an answer to the above ques- 
tion, using non-perturbative methods such as bosoniza- 
tion or the renormalization group (p L G)2iiMiiii£iL2i2iiS. As 
a first step towards a solution of this rather difficult prob- 
lem, it is instructive to study just two coupled spinless 
metallic chain a?^:?: 10 . 

In spite of many years of research, one fundamen- 
tal aspect has only recently been addressed^: the self- 
consistent renormalization of the Fermi surface (FS). The 
band structure of two non-interacting metallic chains 
consists of a bonding and an antibonding band, leading 
in the generic case to a FS consisting of four points ±k F 
and ±kp, as shown in Fig. ^ If the density n is held 
constant, then Luttinger 's theorem implies that the sum 



kp 




FIG. 1: Energy bands of two metallic chains with bare dis- 
persions (antibonding band) and eJL (bonding band), /i 
is the chemical potential. We consider the regime where the 
FS without interactions consists of four points ±kp and ±kp. 
In the text kp and kp always refer to the true FS of the 
interacting system, defined via Eq. Q. 



nn is not affected by interactions; however, 
the difference A = k F — k F , which is proportional to t± 
in the absence of interactions, can be strongly renormal- 
ized. The FS renormalization in two coupled metallic 
chains has been studied previously using RG methods in 
Refsi£*2iiS. However, the calculation of the FS in these 
works is either not self-consistenljSiS in the sense discussed 
in the classic book by Nozieres 11 , or the self-consistency 
has not been properly implemented within the framework 
of the RG™. The self-consistent definition of the FS is 
follows naturally from the fact that the shape of the FS is 
a characteristic property of a zero temperature RG fixed 
p j lr | j i2 1 i3 4 i4 i i5 -yy e em ph as i ze that in our approach the 

concept of a flowing FS depending on the RG flow pa- 
rameter^^ never appears: we construct the true FS as a 
RG fixed point, which by definition does not flow. 

In this work we shall use our general methodiSiiiiiii^ 
of obtaining the FS as a RG fixed point to derive the 
self-consistency conditions for the true Fermi momenta 
k F and k b F of two coupled spinless chains. In the param- 
eter regime where the coupled chain system is a stable 
Luttinger liquid, these conditions can be cast into a sim- 
ple transcendental equation for the distance A = k F — k F 
between the Fermi momenta, which for small t± and for 
weak interactions takes the form 



A = Aill + ^XAo/A) 2 ]" 



(1) 



Here Ai is the value of k F — k F within the Hartree-Fock 
approximation, go is the bare value of a suitably defined 
(see below) dimensionless coupling describing interchain 
backscattering, and Ao 3> A is an ultraviolet cutoff. We 
emphasize that Eq. depends only logarithmically on 
Ao, so that the solution A is not very sensitive to the 
precise numerical value of Ao. Defining x = Ai/A and 
7 = <7q ln(Ao/Ai) 2 , we may rewrite Eq. Q in the form 
x = l+7+.9o hi a; 2 , from which it is easy to sec graphically 
that for 7> 1 even weak interchain backscattering leads 
to a strong reduction of the distance between the two 
Fermi momenta, although for finite 7 they never merge. 
A numerical solution of Eq. Q is shown in Fig. 
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FIG. 2: Distance A = kp — kp between the Fermi momenta 
in the Luttinger liquid phase of two spinless chains coupled 
by weak interchain hopping t±, see Eq. Q. Here Ai is the 
value of k F — k F within the Hartree-Fock approximation, go is 
the bare dimensionless interchain backscattering interaction, 
and 7 = gl ln(A /Ai) 2 . 

Before we discuss of derivation of Eq. (JTJ , let us com- 
ment on its physical significance. First of all, Eq. (JTJ 
predicts a strong "attraction" between the two Fermi 
momenta. It is tempting to extrapolate this result to 
an infinite array of coupled chains: then the FS forms 
a continuum consisting of two corrugated sheets. Our 
two-chain result suggests that the curvature of the sheets 
should decrease due to interactions. This in turn en- 
forces the nesting singularities and hence might stabilize 
a Luttinger liquid state even in higher dimensions. The 
fact that interchain back-scattering reduces the value of 
kp — kp has already been noticed by Fabrizio^. However, 
he neither studied this phenomenon quantitatively, nor 
did he try to calculate the true FS self-consistently. 

The renormalization of k b F — k F predicted by Eq. 
behaves rather different from the estimate based on per- 
turbation theory around the exact solution of the iso- 
lated chains with only forward scatteringLiSii!. Such an 
estimate leads to the prediction that the renormalized 
interchain hopping is t* ± oc t yt ^ / v p K^ 1 ^ , where 77 is 
the anomalous dimension of the Luttinger liquid and vf 
the Fermi velocity for t± = 0. The fact that without in- 
teractions kp — kp oc t± suggests that in the interacting 
system kp — k F oc t j_ oc t±[l — 77 In |i>,fAo/£_l|] for small 
t]. A comparison with Eq. Q shows that this is not a 
valid procedure to estimate the location of the true FS. 

II. EXACT RG EQUATION FOR THE FERMI 
SURFACE 

In order to derive Eq. (JJJ we need to keep track of 
the momentum- and frequency dependence of the four- 
point vertex. To do this we use here the functional RG 
method 12 i 13 . 14 i 18 . 19 i 20 . 21 i 22 i 23 . We have derived Eq. (JTJ 
using both the field theoretical 13 ' 14,21 and the Wilso- 
nian formulationiS^ of the functional RG. In order to 



make contact with previous work^, let us outline here 
the derivation within the Wilsonian approach. We shall 
briefly comment on the corresponding calculation within 
the field theoretical RG in the concluding paragraph. 
Starting point in the Wilsonian approach is the exact 
hierarchy of RG flow equations for the irreducible self- 
energy Sa(ct, K) and the one-particle irreducible 2n- 
point vertices T^ n \{(Ji,Ki\) describing the change of 
these quantities as we vary the infrared cutoff A. Here 
K = (fey , iu>) denotes the momentum ku along the chain 
and the frequency iui. The pseudospin a = a,b labels 
the two bands. We use the Matsubara formalism at zero 
temperature and the same sign conventions for the ver- 
tices as in Ref. (23}. For convenience we introduce the 
cutoff A in momentum space, so that for a given A all 
states with |fc|| ± k F \ > A are integrated out. 

To take the scaling properties close to the Luttinger 
liquid fixed point into account, we introduce rescaled mo- 
menta q and frequencies ie by setting k\\ — ak F + Aq and 
iuj = vpAie, where a — sgnfc|| = ±1 labels the right/left 
half-axis. We define the rescaled two-point vertex by 

T?\cr,a,Q) = -?L[X A (a,K) - E(a,a*#,t0)] , (2) 

which is a function of the rescaled energy-momentum 
Q = (q, ie) and the logarithmic flow parameter / = 
ln(Ao/A). Here Zf is the wave-function renormalization 
factor, and the counter-term £(<x, ak F , iO) is the exact 
self-energy at the exact k F . Once we know the counter- 
term in Eq. J2Jl, we can calculate the true FS for a given 
chemical potential /i from the defining equation^ 

e%*+E(<j,ak F> iO)=iJ,, (3) 

where ejL is the energy dispersion of band a = a,b. 
In higher dimensions the corresponding counter-term 
S(kj?,i0) and the momenta ki? on the FS can be cal- 
culated within self-consistent perturbation theory in the 
weak coupling regime 24 . On the other hand, in one di- 
mension it is necessary to invoke the RG to properly treat 
the infrared divergencies. A formally exact equation for 
the counter-term can be obtained from the requirement 
that for I — > 00 the relevant part 

rf EElfV,a,0) (4) 

of the two-point vertex (J2J flows into a RG fixed points. 
Defining the rescaled four-point vertex via 

f\ 4 \a[a[Q[,<j' 2 a' 2 Q' 2 ;cr2Ct2Q2,cr 1 a 1 Qi) = (2Trv F y 1 

i=l,2 

~ (2) 

the exact RG flow equation for T) '(cr,a,Q) reads 

9 ; f ; (2) = (1 - rf[ - Q ■ d Q )r[ 2) + tf\aaQ) , (6) 
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where rff = —di In Zf is the flowing anomalous dimen- 
sion, Q ■ 8q — qd q + ed e , and the inhomogeneity is 

tf\aaQ) = -Y.j d 4%G l {*i a >Qi) 

a' 'a' 

xf j 4) (aaQ, a'a'Q'; a'a'Q', aaQ) . (7) 

Here for sharp cutoff and linear dispersion 

G^aaQ) = -6(l-\q\)/[Zf{ie-av^q)-tf\aaQ)] , (8) 

with the dimensionless Fermi velocities Vq — v F /vp, 
where v F is the bare Fermi velocity associated with band 
a. Because the coupling rf is relevant with scaling di- 
mension +1, its initial value Tq must to be fine tuned in 
order to approach a finite limit for I — > oo. This leads to 
the condition^ 

/■oo 

rl = - dle- l +H^tl 2 \a,a,0), (9) 
Jo 

which relates the inital value to the RG flow on the 
entire RG trajectory. Assuming that Ao is sufficiently 
large so that £a can be neglected and Zq m 1 , we obtain 
from Eq. J5J for the counter-term 

T,{a,ak F ,iO) = ~v F A n r% . (10) 

At constant density n — (k F + k F )/ir we then find 

k b F - k a F = A = A + 2Ao(«g + wS) _1 (rg - r Q ) , (11) 

where Ao is the value of A;^ — k F at the same density but 
without interactions. 

III. WEAK COUPLING TRUNCATION FOR 
THE EFFECTIVE INTERACTION AND THE 
SELF-ENERGY 

The formally exact Eqs. ijfil lllll are now our starting 
point to calculate A self-consistently. Note that the con- 
cept of a "flowing FS" never appears in our approach; 
by definition, k F and k b F are the true Fermi momenta of 
the interacting system, to be determined self-consistently 
from Eqs. (JHJ) and l|ll(l . To make further progress, we 
need an approximate expression for the effective interac- 
tion f in Eq. J7J). To simplify the analysis, we neglect 
vertices describing intrachain umklapp scattering as well 
as chiral vertices, involving particles moving in the same 
direction. The latter are expected to give finite renor- 
malizations of the Fermi velocities, which we ignore here. 
Totally we should then keep track of the RG flow of the 
following five vertex functions, 

jf (Q'i, Q' 2 \ Q a , Qx) = f ; (4) (&Qi, &Q' 2 ; bQ 2 , 6Qi),(12a) 
friQ'v Q' 2 ; Q2,Qx) = f | 4) (oQi, aQ' 2 ; aQ 2 , oQi),(12b) 
ftiQ'^Q'v Qa, Qi) = f | 4) (aQ;, bQ' 2 ; bQ 2 , aQi), (12c) 
MQ'i, Q 2 , Q2, Qi) = fl 4) (bQ' v bQ' 2 ; aQ 2 , oQi),(12d) 
gi(Q'i,Q' 2 ; Q2, Qi) = f l 4) {bQ[, aQ' 2 ; bQ 2 , aQi), (12e) 



where on the right-hand side the direction labels are 
(a[, a' 2 ; <*2, oil) = (+> — i — >+)• Here / ; <Tcr describes for- 
ward scattering of two fermions associated with bands a 
and a', ui describes interband umklapp scattering (also 
called pair tunneling), and gi corresponds to interband 
backscattering. The functional RG equations for these 
five functions follow from the general flow equations given 
in Ref; 1 -- in a straightforward way. If we set Qi = Q\ = 
and denote the corresponding coupling constants by fj 717 , 
ui and gi (without a tilde), we obtain the one-loop RG 
flow equations 

dm = -uifi , digi = gifi , difi = -Auf + B t gf , (13) 

where 

fi = frn + jf m - vt b /v , (w) 

with A = 2/(vqVq) + 2/vq and the average dimensionless 
velocity vq = (vq + Vq)/2. The function 

Bi = 2C , ,/(«5«J) + 2C,/wg (15) 

depends on the rescaled Fermi point difference A; = 
Ae l /A via 

Q = 6(e i -l-2|A i |)/(l + |A i |) , (16) 
Oi = 8(e'-l-2|A ; |)i^ _ * (17) 

The RG equations (|13|l are equivalent to those derived 
by Fabrizio^ at one-loop order. Note that the RG flow of 
ui and gi couples to the forward scattering interactions 
only via the combination /; defined above. The flow of 

ir' is 

tW, 66 - -uf + Qgf = v b d t fr , (18) 
%dift = «? - C l9 f . (19) 

As discussed by Fabrizio^, the above RG equations pre- 
dict a finite regime of initial values f$ a , Uo, and go where 
the system is a stable Luttinger liquid, characterized by 
finite forward scattering couplings // Tcr and vanishing ui 
and gi for I — > oo. 

We now calculate the true FS in the Luttinger liquid 
regime using our general Eqs. ijBT lllfl . Assuming that 
fi* 7 , ui, and gi are small compared with unity, we may 
adopt the same strategy as in Refsi22i2£: we expand the 
right-hand side of the flow equations for the momentum- 
and frequency dependent vertices defined in Eqs. (|12al - 
I12e|) to second order in fj 7 ' 7 , ui , and gi . The flow of the 
above vertex functions is then given by the following set 
of equations, 

D bb f bb (Q' 1 ,Q' 2 ;Q 2 ,Q 1 ) = 

{ft b ) 2 [xT{Qi-Q' 2 )~xT{Qi + Q 2 )\ 

-uh? a (Qi + Qi) + 9iX? a (Qi ~Q' 2 + 2A Z ) , (20a) 
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where for a,(3,j G {a,b} the integrals /"f (Q) are 



Dt a fr(Q' 1 ,Q' 2 ;Q2,Qi) = 

(fn 2 ixr(Qi-Q2)-xr(Qi + Q2)} 

-uf X ?(Qi + Qa) + & 2 Af(Qi - Qa - 2A . ( 2 0b) 



Dfft h (Q[,Q' 2 ;Q 2 ,Qi) = 

(f? b ) 2 lxf(Qi~Q2)-xf(Qi 

+ufx1 a (Qi-Q'2)~9?x b i a (Qi- 



Dfu l (Q' 1 ,Q' 2 ;Q 2 ,Q 1 ) 



-Qa)] 

Q 2 - 2A ; ) , (20c) 



/r^[xf(Qi 



Q 2 )+xr(Qi-Q 2 )] 

o 2 ) - frnxriQi - 



2), (20d) 



A ab 3;(Q'i,Q 2 ;Q2,Qi) = 

-/r h <7ax? b (Qi + Q2)+xf a (gi 



2)] 



+fimxi(Qi - Q' 2 ) + frmxt a (Qx - Q 2 ) • (20e) 

Here Q + A; is a short notation for (g + A/, ie), the gen- 
eralized susceptibilities are given by 



xf(Q) 



&(e l 



1 



M) 



2t>o + «o M - * £ sgng 

e(e ; - 1 - M) 



2«o 



lesgnq 



(21) 



where a,/3G {a, 6}, and 



0°' 



9; + + »?f + £ W< ' + " ^0< ) (22) 



are comoving derivatives. The linear partial differential 
equations H20all20ej) can easily be solved exactly 2 ^, so 
that we can obtain an explicit expression for the flow- 
ing momentum- and frequency dependent effective inter- 
action T^ 4 \aaQ, a'a'Q'; a'a'Q' , aaQ). Substituting this 
expression into Eq. (JJJ and performing the Q'-integration 
we find 



t{ 2) (a,+,Q) = f? a + n b - di[fr + f? b \ 
-2(/D 2 C(0) - 2(ft) 2 iti(Q) 

-uf[i%i(Q)+iftm 

-gf {C (q-2A h ie) + 1$ (q - 2 A, , ie)} , (23) 



V Q n =±l 



nQ(e l -1- \q + n\) 



x |6(ns g +„) In 
+@(-ns q+n )ln 







+ ^0 


-< 


9 


+ rt 1 — me 






«0 + 


"ok 


+ 


n\ — ine 






+ v y 






+ rt | — me 










+ 


77 1 — 7776 



(24) 



with s g+ „ = sgn(q + n). The function T ; (&, +, Q) can 
be obtained from Eq. (|2*3)) by replacing a «-> 6 and A; — * 
—A; on the right-hand side. Once we know the function 

Ij (er, +,<5), we may calculate the flowing anomalous 
dimension from 2 ^ 

r / r = -5ff ) ( ( x,+,Q)/9(7 e )| Q= o. (25) 
We can reproduce the result by Louis et al~ if we sim- 



ply set f{ 2) (a, 



,Q) « /,«» + /f b instead of Eq. J2U, 



which amounts to replacing the effective interaction 
(aaQ, a'a'Q'; a'a'Q', aaQ) in Eq. Q by its value 
at Q = Q' = 0. This approximation does not consis- 
tently take into account all contributions to rg which are 
quadratic in the bare couplings, because such second or- 
der terms are also generated by the Q'-dependence of the 
four-point vertex in Eq. (0. 

To calculate the FS from Eq. Qllfl. we still need to 
perform the trajectory integral 10. Since we are inter- 
ested in the leading behavior for small t±, we may at 
this point ignore the small difference between the bare 
Fermi velocities, setting Vq = Vq = vq = 1. To sec- 
ond order in the couplings, we may also set 77J 7 = in 
Eq. 10. With these approximations the integral obtained 
by substituting Eq. (|2*3)) into Eq. © can be performed 
analytically to second order in the bare couplings. For 
A = A/A < 1 the result is rf} = -f$ a - / ab - 5 2 A In | A| 
and r b = -f$ b - f ba +^Aln|A| to leading logarith- 
mic order. Substituting this into Eq. i|f f [I we arrive at 
Eq. Q, with Ax = A - A (/ bb - /o Qa ), where we have 
used fS b = / bQ . 

Finally, let us point out that we have also derived 
Eq. 0} within the field-theoretical RG, where we calcu- 
late the renormalized vertex functions directly from per- 
turbation theory. Setting again Vq — Vq — 1 we obtain for 
the difference between the self-energies = Y,(a, kp,i0) 
at the true FS, up to two-loop order, 



-2f 1 



v F Ao[f b R b ~f a R a -g 2 R A\nA 2 } , (26) 



where /j^ CT and g R are the renormalized couplings. Since 
up to this order we do not distinguish between bare and 
renormalized couplings, the combination of this result 
with the definition J3J immediately leads to Eq. (JTJ. 



5 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have presented a fully self-consistent 
calculation of the true FS of two coupled metallic chains 
in the regime where the Luttinger liquid fixed point is 
stable. Our final result for the renormalized FS is given 
by the solution of the transcendental equation Q shown 
in Fig. which is perhaps the simplest example for an 
explicit solution of the self-consistency problem for the 
true FS discussed by Nozieresii. We find that for small 
interchain hopping t± even weak interchain backscatter- 
ing can lead to a strong reduction of the distance be- 
tween the Fermi points associated with the bonding and 



the antibonding band. It is tempting to speculate that in 
higher dimensions a similar smoothing effect of a corru- 
gated FS stabilizes a strongly correlated non- Fermi liquid 
state. Our method is general and can be used to calculate 
the true FS in higher dimensions within the framework 
of the RG; of course, in this case the resulting integral 
equations can only be solved numerically. 

This work was completed while P.K. was visiting the 
International Center of Condensed Matter Physics (IC- 
CMP) at the University of Brasilia, Brazil. He would like 
to thank the ICCMP for its hospitality and for financial 
support. 
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